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*^ , We calculate the total energy of Taub's 1951 exact solution for a Bianchi type- 

■ IX geometry using several different energy-localization procedures, including the 

prescriptions of Einstein, Papapetrou, Landau-Lifshitz and M0ller. We compare 
these results to those for other anisotropic geometries, and comment on their 
relationship to Rosen's conjecture about the total energy of the universe. 



3 



(N 

(N 
> 

in 

O 

\o 

O 



- 1—1 

X 



Keywords: Taub solution; Mixmaster universe; energy-momentum complex. 
PACS numbers: 04.20.Cv, 04.20.Jb, 98.80.Jk 



1. Introduction 



A longstanding question in general relativity is whether gravitational energy 
can be localized and, if so, what measure to use. Einstein proposed the first 
\Q \ energy-momentum complex in an attempt to define the local distribution of 

energy and momentum [1]. Einstein's measure is non-tensorial and asymmetric 
O i in its indices. These issues have inspired a number of alternative measures, 

including prescriptions by M0ller [2], Tolman [3], Papapetrou [4], Landau and 
£h ■ Lifshitz [5], Weinberg [6] and so forth. By in large, these formulations depend on 

the coordinate-system used. For example, Einstein's prescription favors quasi- 
Cartesian coordinates, defined as those coordinates x 1 for which the metric 
components converge sufficiently rapidly toward the constant components 
of the diagonal Minkowski metric r/^. Some of these complexes have other 
restrictions. 

Such limitations have led some researchers to discount the importance of 
localized measures. However, in 1990 Bondi proved that energy must be localiz- 
able, at least in principle [7]. The same year, Virbhadra demonstrated that these 
complexes, though non-tensors, yield credible, consistent results when applied 
to the same spacetime geometries [8]. In 1993, Rosen and Virbhadra employed 
the Einstein prescription and found reasonable results for the energy and mo- 
mentum of cylindrical gravitational waves [9]. Later, Virbhadra derived similar 
results using Tolman's procedure and the measure proposed by Landau and Lif- 
shitz [10]. Then in 1996, Aguirregabiria, Chamorro and Virbhadra found that 
for any metric of the Kerr-Schild class (including the Schwarzschild, Reissner- 
Nordstrom, Kerr and Vaidya metrics, for example.), various energy-momentum 
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prescriptions yielded reasonable and consistent results [11]. Of these, Virbhadra 
found that the Einstein prescription offered the greatest consistency [12]. These 
findings have stimulated a number of investigations attempting to extend such 
methods to other classes of metrics [13]. 

Much work has been done, for example, in using various measures to deter- 
mine the energy distributions of assorted generalizations of black holes. Chamorro 
and Virbhardra [14] as well as Xulu [15] have examined the energy distribution of 
charged black holes with a dilaton field. Radinschi and Yang [16], Vagenas [17], 
Gad [18] and Xulu [19] have each looked at the energy distribution of stringy 
black holes. Vagenas has also examined black holes in 2+1 dimensions [20, 21] 
and applied M0ller's prescription to the dyadosphcrc of a Reissner-Nordstrom 
black hole [22] . The energy distribution of black plane solutions has also been 
investigated [23]. Other researchers have looked at the energy distribution for 
the case of teleparallel gravity [24], for a stationary beam of light [25, 26], for a 
radiating charged particle [27] and for a host of other metrics. 

Cosmological solutions are another area of considerable interest, stimulated 
by a conjecture made by Rosen in 1994 that the total energy of the universe, 
presuming it has a closed geometry, is zero [28] . Applying the Einstein energy- 
momentum complex to a closed Friedmann-Robertson- Walker (FRW) cosmol- 
ogy, Rosen determined its total energy to be indeed zero. Johri, Kalligas, Singh 
and Everitt later reached similar conclusions by applying the prescription of 
Landau and Lifshitz [29]. 

One possible implication of Rosen's conjecture is a hypothesis proposed by 
Cooperstock that the energy distribution is non-zero only in regions where the 
stress-energy tensor T M „ is non-zero [30] . This implies that gravitational waves 
traveling through a vacuum would carry only information content, not energy, 
and thus would be fundamentally different from other forms of radiation. 

Researchers have sought to extend these results to more general cosmologies. 
The isotropic FRW model is far from the most general cosmological solution, 
and therefore would not, for example, represent a chaotic very early stage of 
the universe. The simplest anisotropic model is the Kasner solution, possess- 
ing a Bianchi type I geometry. Several researchers have investigated the total 
energy of this model using a variety of measures, and have found the result 
to be identically zero for any closed region, consistent with Rosen's conjecture 
[31, 32, 33, 34, 35]. In addition to the standard vacuum case, Salti looked at 
the viscous Kasner universe in teleparallel gravity [36] using the Einstein and 
Landau-Lifshitz formulations and found the total energy and momentum to be 
identically zero. Generalizing these findings to another anisotropic cosmology, 
Radinschi examined the energy-momentum of the Bianchi type VIO universe 
using a number of different measures, finding that all the components vanish 
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[37]. 



2. The Taub Cosmological Solution 

Amongst homogenous cosmologies, the metric that yields the most complex 
anisotropic behavior possesses a Bianchi type IX spatial geometry. The one- 
forms associated with the spatial part of this metric possess the Lie algebra 
associated with one of the normal subgroups of the group of motion of a three- 
dimensional sphere. 

As Belinski, Khalatnikov and Lifshitz demonstrated in 1970, in evolving the 
vacuum case backward in time toward its initial singularity, the behavior of 
its three scale factors resembles a sequence of transitions between Kasner-likc 
epochs [38]. In using Hamiltonian methods to analyze this behavior, Misner 
dubbed it the "Mixmastcr universe," because of its special mixing properties 
[39]. Barrow later demonstrated that the vacuum Bianchi type IX solution 
exhibits a form of deterministic chaos [40] . 

If one constrains two of the three scale factors to be equal, one obtains an 
exact, non-chaotic solution, first identified by Taub in 1951. The Taub metric 
can be expressed in a holonomic coordinate system as: 



ds 2 = dt 2 - R(t) dx 2 - (R(t) sin 2 (x) + S(t) cos 2 (x)) dy 2 

- 2 S(t) cos(x) dy dz - S(t) dz 2 (1) 

with x, y and z corresponding to the Euler angles of three-dimensional ro- 
tation, and having the ranges: 

< x < tt; < y < 4tt; < z < 2ir (2) 

As Taub demonstrated, to satisfy the vacuum Einstein equations the scale 
factors are constrained to be: 



_ k cosh(fc t + a) 
R{t) 4 cosh 2 (^) (3) 

SW - cosh( fc \ + .) (4) 



where a, (3 and k are constants and with: 



7 (t) = R(tfS(t) (5) 
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One reason the Taub solution is of interest is that in 1968, Misner and Taub 
[41] proved that it represents the spacetime region bounded by the horizon of 
the NUT solution found by Newman, Unti and Tamburino [42]. Joined together, 
these constitute the Taub-NUT spacetime. 



3. Finding the Energy Using Einstein's Energy-Momentum Complex 

It is interesting to calculate the total energy of the Taub cosmology using 
several different measures in order to make a comparison of the results. We 
begin with the Einstein prescription, for which the energy-momentum complex 
is defined as: 

e = i(b^ (6) 

with the superpotentials H* 1 given by: 

H M = Jhr^ hg{g kn g l m _ gln g km )] m (?) 

which possesses the antisymmetric property that: 

H t M = -Hi lk (8) 

Substituting into (7) the metric components of the Taub solution (1-5) we 
find that the only non-zero components of H kl are: 

H 01 = -H la = cos(x) sech 2 (^±^) (9) 

o Z 

Therefore, from (6) the only non-zero components of 9i k are: 

^o 1 = cos(x) sech 2 (^) tanh(^) (11) 

In general, the energy-momentum components can be found by the volume 
integral: 

P t = J J J 9 i °dx 1 dx 2 dx 3 (12) 
Thus, if we substitute the coordinate ranges (2) we can write the total energy 

as: 



4 



P Q = / 9 °dx dy dz (13) 

Jo Jo Jo 

Evaluating this integral for (10), we obtain an expression for the total energy 
of the Taub solution: 

* = Po = ^sech2 ( ^) (u) 



4. Following Papapetrou's Procedure 

Unlike Einstein's prescription, the energy-momentum complex of Papapetrou 
has the advantage of being symmetric in its indices. Hence it allows local con- 
servation laws to be well-defined. 

The Papapetrou energy-momentum complex is defined as: 

W k = —N**** ab (15) 

Pi = J J J n t0 dx 1 dx 2 dx 3 (16) 
where the functions N tkab are given by: 

N lkab = [g ik r] ab - g ia r] kb + g ab rf k - g kb rf a \ (17) 

and the g ab terms represent the components of a Minkowski metric of sig- 
nature —2. 

Inserting the metric components of the Taub solution (1-5) into expression 
(17), we obtain the following relevant non-zero components of N lkab : 

N oou = N uoo = _ sin(x) [4 + k 2 sech 2 { kt+£y (lg) 

jyOlOl = jylOlO = gin(x) [4 + k 2 scch 2 { '^+l )] (19) 

Substituting into (15), we find that the only non- vanishing component of the 
energy density is: 

Q00 = -I- sin(x) [4 + k 2 sech 2 (^^)] (20) 
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To obtain the total energy, we integrate the energy density f2 00 over the 
appropriate coordinate ranges (2): 



Pq= / tt°°dx dy dz (21) 

Jo Jo Jo 

This yields: 



E = P = 7T [4 + k 2 sech 2 ( fc ^ + /? )] (22) 



5. Using the Technique of Landau and Lifshitz 

Like Papapetrou's procedure, the Landau-Lifshitz energy-momentum com- 
plex is symmetric in its indices, permitting conservation laws, such as conserva- 
tion of angular momentum, to be readily defined. 

We define the complex by means of the functions: 



L" = i^,r (23) 
where the L t0 represent the energy and momentum densities and: 



S W = _ g ^gkl _ gikgjl} (24) 

Substituting into this expression the metric components of the Taub solution 
(1-5) we find the only relevant non-zero component to be: 

S 0011 = -^sm 2 ( x) S ech 2 ( ^) (25) 
From (23), we determine the energy density L 00 to be: 

L ™ = -^ cos(2x) sech 2 (^±^) (26) 
We now integrate this expression over the full range of coordinates: 

Po= L m dx dy dz (27) 

Jo Jo Jo 

Evaluating this integral, we find the total energy to be: 

E = P a = (28) 
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So far, we have seen three different results for the total energy, from the 
Einstein, Papapetrou and Landau-Lifshitz complexes respectively. Perhaps 
this discrepancy stems from the coordinate system used, which is not quasi- 
Cartesian. To obtain coordinate-independent results, we now turn to Miller's 
energy-momentum complex, specially designed to avoid such issues. 



6. Using M0ller's Method 

M0ller constructed his complex to transform spatially like a tensor and hence 
be as free as possible of the choice of coordinate system. Hence, unlike other 
complexes, it does not favor quasi-Cartesian coordinates. Given that the nat- 
ural coordinates to define the Taub solution arc not quasi-Cartesian, Miller's 
procedure seems an ideal method to evaluate its gravitational energy. 

The quantity M0ller developed is defined as: 



ixt (29) 

where: 



xT = V~g [g iP , q g iq , P } g kq g lp (30) 

Inserting the Taub metric components (1-5) we find that for all values of k 
and I: 



Xt l = (31) 



Therefore, for all values of k: 



"0 

Integrating, we obtain the total energy: 



7T />47T f>2lT 

"0 



" k - (32) 



E= / / E£ dx dy dz = (33) 
Jo Jo Jo 

Thus, M0llcr's energy-momentum complex agrees with the prescription of 
Landau-Lifshitz that the total energy of the vacuum Taub cosmology is zero. 
We also find that all the momentum components are zero. 



7. Conclusion 
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We have investigated the total gravitational energy of the vacuum-filled Taub 
cosmology using four different energy-momentum complexes: those of Einstein, 
Papapetrou, Landau-Lifshitz and M0ller. Of these, the latter two yielded a value 
of zero for the total energy of the universe, agreeing with Rosen's conjecture that 
the total energy of any closed cosmology is zero. The Einstein and Papapetrou 
complexes, in contrast, yielded non-zero quantities that evolve over time. In 
our opinion, these two results are suspect because the natural coordinates for 
the Taub metric, though holonomic, are not quasi-Cartesian, and thus do not 
fulfill the ideal conditions for these complexes. 

We note that our results from the M0ller and Landau-Lifshitz complexes 
agree with earlier work on anisotropic cosmologies, including evaluations of the 
total energy of the Kasner (Bianchi type I) cosmology for several different cases. 
We conjecture that the total energy of all vacuum-filled anisotropic cosmologies 
is zero. Further calculations, including other Bianchi types, would be needed to 
confirm this hypothesis. 
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